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, Abstract We consider a finite volume scheme for the two-dimensional in- 

^— i ■ compressible Navier-Stokes equations. We use a triangular mesh. The un- 

knowns for the velocity and pressure are respectively piecewise constant and 
affine. We use a projection method to deal with the incompressibility con- 
straint. The stability of the scheme has been proven in [15]. We infer from it 
its convergence. 
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1 Introduction 



> 

oo 
o , 

■ ' | We consider the flow of an incompressible fluid in a open bounded polyhedral 

set Q C M 2 during the time interval [0, T}. The velocity field u : Q x [0, T] -> 
M 2 and the pressure field p : Q x [0, T] — > R satisfy the Navier-Stokes equa- 
f^*) . tions 



u t --^-Au+(u-V)u + Vp=f, (1) 
Re 

div u = , (2) 

with the boundary and initial condition 

u|afi = 0, u| 4= o = u . 

The terms Au and (u • V)u are associated with the physical phenomena of 
diffusion and convection, respectively. The Reynolds number Re measures the 
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influence of convection in the flow. For equations {J}-©, finite clement and 
finite difference methods are well known and mathematical studies are avail- 
able (see [9] for example). For finite volume schemes, numerous computations 
have been conducted ([T3] and [T] for example). However, few mathematical 
results are available in this case. Let us cite Eymard and Herbin [5] and 
Eymard, Latche and Herbin [6]. In order to deal with the incompress- 
ibility constraint, these works use a penalization method. Another way is 
to use the projection methods which have been introduced by Chorin [3] 
and Temam [15] . This is the case in Faure [8] where the mesh is made of 
squares. In Zimmermann [14] the mesh is made of triangles, which allows 
more complex geometries. In the present paper the mesh is also made of tri- 
angles, but we consider a different discretisation for the pressure. It leads to 
a linear system with a better-conditioned matrix. The layout of the article is 
the following. We first introduce (section |2~T]) some notations and hypotheses 
on the mesh. We define (section 12 . 2[) the spaces we use to approximate the 
velocity and pressure. We define also (section 12. 3[) the operators we use to 
approximate the differential operators in (|TJ) — d2I) - By combining this with a 
projection method, we build the scheme in section [3J In order to provide a 
mathematical analysis, we state in section [4] that the differential operators in 
((T])-([2]) and their discrete counterparts share similar properties. In particular, 
the discrete operators for the gradient and the divergence are adjoint. We 
then prove in section [5] the convergence of the scheme. 

We conclude with some notations. We denote by \i the characteristic func- 
tion of an interval I cK. We denote by = C^°(J7) the set of the functions 
with a compact support in fi. The spaces (L 2 ,|.|) and (L°°, ||.||oo) are the 
usual Lebesgue spaces and we set Lq = {q e L 2 ; J„ q dx = 0}. Their vectorial 
counterparts are (L 2 , |.|) and (L°°, y^) with L 2 = (L 2 ) 2 and L°° = (L°°). 
For k e N*, (H k 1 \\ ■ \\ k ) is the usual Sobolev space. Its vectorial counterpart 
is (H fc , ||.|| fe ) with H k = (H k ) 2 . For k = 1, the functions of H 1 with a null 
trace on the boundary form the space Hg. Also, we set Vu = (Viti, V«2) T 
if u = (ui,U2) G H 1 . If X C L 2 is a Banach space, we define C(0,T;X) 
(resp. L 2 (0,T;X)) as the set of the applications g : [0, T] — > X such that 
t — > |g(t)| is continous (resp. square integrable). The norm ||.||c(o,T;X) is 
defined by ||g||e(o,T;X) = su Pse[o,T] |s( s )|- Finally in all calculations, C is a 
generic positive constant, depending only on J?, Uo and f. 

2 Discrete setting 

First, we introduce the spaces and operators needed to build the mesh. 
2.1 The mesh 

Let Th be a triangular mesh of Q: fl = UKeT h K. For each triangle K £ Th, 
we denote by \K\ its area and £k the set of his edges. If a e £k, xik,<j is the 
unit vector normal to a pointing outwards of K . 

The set of edges of the mesh is Eh = U k eT h £ k ■ The length of an edge a € £h 
is \a\. The set of edges inside Q (resp. on the boundary) is £ h nt (resp. £| xt ): 
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£ h = £™ U ■ If c € £h ' K a and L a are the triangles sharing a as an 
edge. If a £ £^ x *, only the triangle K a inside Q is defined. 
We denote by the circumcenter of a triangle K . We assume that the 
measure of all interior angles of the triangles of the mesh are below so 
that xk G K. If a £ £™* (resp. a £ £% xt ) we set d a — rffxjc^Xi^) (resp. 
d a = d(x a , x^)). We define for all edge c G tv = j^. The maximum 
circumradius of the triangles of the mesh is h. We assume that there exists 
C > such that 

Va££ h , d(x Ki7 ,(T) > C\a\ and \a\ > Ch. 

It implies that there exists a constant C > such that for all edge a £ £h 

r a >C (3) 

and for all triangles K £ Th we have (with a £ £k and Kk.o the matching 
altitude) 

\K\ = I M h K , a > \ \a\ d(x K ,x CT ) > Ch 2 . (4) 



2.2 The discrete spaces 
We first define 

Pq = {q £ L 2 ; M K £ %, q\ K is a constant} , P = (P ) 2 . 

For the sake of concision, we set for all qh £ Po (resp. Vft £ Po) and all 
triangle K £ T h : q K = qh\x (resp. w K = vh\k). Although P (f. H 1 , we 
define the discrete equivalent of a H 1 norm as follow. For all v/i £ Pq we set 



K|k= £ r ff |v L „-v K J 2 + £ ^|v*J 2 • (5) 




We have [7] a discrete Poincare inequality for Po: there exists C > such 
that for all v/ t £ Po 

|v fc |<C|K|U. (6) 
From the norm \\.\\h we deduce a dual norm. For all v/, £ Po we set 

huU,H= sup ^1^1. (7) 

For all u/j G P and Vft G P we have (u/^v/,) < Hu/jH-i^ \\v h \\ h . We 
define the projection operator 77p : L 2 — > Po as follows. For all w G L 2 , 
i7p w G Po is given by 

Vif G 7ft , (i7 Po w)|x = 7-^7 J w(x) dx. (8) 
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We easily check that for all w G L 2 and v/ t G Po we have (JJp w,v^) = 
(w, v/j). We deduce from this that 7Ip is stable for the L 2 norm. We define 
also the operator 77p : H 2 — > Po- For all v 6 H 2 , -/7p v 6 Po is given by 

VKeT h , n Po v\ K = v{x K ). 

According to the Sobolev embedding theorem, v € H 2 is a.e. equal to a 
continuous function. Therefore the definition above makes sense. One checks 
[14] that there exists C > such that 

|v-i7p v| <Ch\\v\\ 2 (9) 

for all v G H 1 . We introduce also the finite element spaces 

Pf = {v G L 2 ; VK G T h , v\ K is affine} , 

Pl lc = {vh e Pf ; V, e S h nt , v h \ Ka {^) = v h \ L „(* a ) , P? c = (PH 2 - 

If qh G -P" c , we have usually Vg^ ^ L 2 . Therefore we define the operator 
V/j : P™ 5 -* P by setting for all qt G Po an d all triangle K e Th 

V h q h \K = t4| / V * dx - ( 10 ) 
We define the projection operator Ptp^c. For all q £ if 1 , Ilp^^q is given by 

Vff£4, / (flpf?) ^ = f qda. (11) 

We also set TIp^c = (Pp™c) 2 . One checks that there exists C > such that 
|Vg-V fc (17j». ? )|<C/i|| g ||2 > |v-J7p r v| <CA||v||i, (12) 

for all q G H 1 and v G H 1 . 

We also use the Raviart-Thomas spaces 

RTq = {v ft G Pf ; Vct G Vft|if • njf^ is constant, and v h ■ n\ d a = 0} , 
RT = {v ft G RTJJ ; VK G T ft , Vct G S K , v h \ Ka ■ n Ka>(T = v h \ L „ ■ n K „, a }. 

For all v h G RT , K G % and cr G £r- we set (v h ■ n K . a ) a = \ h \ K ■ n Kt<J . 
We define the operator P RTo : H 1 -> RT . For all v G H 1 , 7T RTo v G RT 
is given by 

VKeT h , Vctg£k, (n RTo v ■ n KM ) a = -- / vda. (13) 

One checks [3] that there exists a constant C > such that for all v G H 1 

|v-i7 RTo v| < C/»||v||i. (14) 
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2.3 The discrete operators 

The equations {l}-© use the differential operators gradient, divergence and 
laplacian. Using the spaces of section \2?2\ we now define their discrete coun- 
terparts. The discrete gradient V/, : P[ iC — ► Po is defined by ([TO]) . The 
discrete divergence operator div^ : Po — > -P" c is built so that it is adjoint 
to the operator Vft (proposition [3] below) . We set for all € P and all 
triangle K E Th 

3 \o\ 

e 4"*= (divh v h )(x ff ) = i l ™ i (v Zff - v*J • n KjCT ; 

Vaefr*, (div h v fo )(x ff ) = - J ^ Vgg • n K , a . (15) 

The first discrete laplacian Ah ■ P\ c — > -P" c is given by 

V % eP; c , 4 ft % = div h (VhQh)- 

The second discrete laplacian : P — > Po is the usual operator in finite 
volume schemes [7]. We set for all v/j 6 Po and all triangle K £ T h 

A h y h \ K = r^7 t °( w l« ~ - T^T H T ^ v ^- (16) 

In order to approximate the convection term (u • V)u in (fT]), we define a 
bilinear form : RTo X Po — * Po using the well-known [7] upwind scheme. 
For all u/j e Po, G Po, and all triangle K G Th we set 

Mufc.v/,)^ = r^r ^ |ct| ((u • nj£- )(T )+ v K + (u ■ n K><7 )- v L<r ). 

We have set a + = max(a, 0), a - = min(a, 0) for all net. Lastly, we define 
the trilinear form bh ■ RT x P x P — > K as follows. For all Uh e RT , 
V/, € P , w,j € P , we set 

b h (uh,, Vh, w ft ) = ^ \K\w K ■b h {u h ,v h )\ K . 

KeT h 



3 The scheme 

We have defined in section [5] the discretization in space. We now have to 
define the discretization in time, and treat the incompressibility constraint 
We use a projection method to this end. This kind of method has been 
introduced by Chorin [3] and Temam [T3j. The time interval [0, T] is split 
with a time step k: [0,T] = U n =o[*»' ^n+i] with N e N* ct t n = nk for all 
n E {0, . . . , N}. We start with the initial values 

u° e Po n RTo , u\ e P n RT , y\ e Pr c n Lq . 

For all n £ {l,...,N}, (u£ +1 ,p£ +1 , u£ +1 ) is deduced from (u£,p£,u£) as 
follows. 
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- < +1 G Po is given by 

3< +1 -4< + <- 1 1 ~ +1 
2k Rc AhUh 

+b„(2ug - v^k - fr 1 , (l?) 

- pl +1 G P'[ lc n L\ is the solution of 

^ +1 -^) = ^div,u«+ 1 ! 

- u^ l+1 G P is given by 

K +1 = K +1 - 2 ^ h (vl +1 -P^- (is) 

We have proven in |14j that the scheme is well defined. In particular the term 
b/j(2u^ — u^ _1 ) in (|17| is defined thanks to the following result. 

Proposition 1 For m G {0, . . . , N} we have u™ G RTo- 

Note also that for m G {0, . . . , N} we have divu™ = 0, since u™ G Po- Thus 
the incompressibility condition {2} is fullfilled. 

4 Properties of the discrete operators 

The operators defined in section 12.31 have the following properties [M] . 

Proposition 2 There exists a constant C > such that for all u/j G RTo 
satisfying divu,j = 0, v h G P , w h G P : 

Ibfc^.VftjVfc)! < C|ufe| ||vfe|| ft HwftHft. 

Proposition 3 For allw h G P and % G -P" c .- (vfe,Vft%) = -(<&, divft, v ft ). 

Proposition 4 For all G Po and v h G Po-' — (A^u^jU/,) = ||uft||^ and 
-(A^UftjV/,) < \\uh\\h \\yh\\h- 

If v G H 1 we have |divv| < ||v|ji. The operator divh has a similar property. 
Proposition 5 There exists a constant C > such that for all Vh G Po 

|div/,v h | < C\\v h \\ h . 
Proof. Using a quadrature formula we have 

|div h v,| 2 = E ^ E l(div/ l v,)(x CT )| 2 . 

KGT h cr££ K 
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Let K £ Th- Using definition and ((4]) we have 

Va e £ K n 4 ni > |(div h v, i )(x CT )| 2 < ^ |v Lct - v^| 2 ; 
Vae^nfr', |(div^v,)(x CT )| 2 < ^|v^| 2 . 

Thus: |div fc v h | 2 < CJ2 KeT [j2 <7££Kr]£rt \v La - v K \ 2 + E a e£ K n£^W K \ 2 ) ■ 
Writing the sum over the triangles as a sum over the edges, we get 

\div h v h \ 2 <C I r a \v L<r -v K \ 2 + ]T r a \v K A <C\\v h \\l m 

Proposition 6 Ifu h S Po andv h £ Po we We (A/jii/j, v;,.) = (u ft , A^v^). 
PROOF. Using definition (|2.3|) one checks that 

(AftUft,, v fe ) = ^ (v L<T - v K J ■ (u i<T - urtJ - ^ r CT v K<T • \i Kcr 

= (A h v h ,u h ). m 

Proposition 7 There exists C > swc/i that for all v e H 2 satisfying 
Vv ■ n| 9fi = 

||77 Po (Av)- A h (fl- Po v)||_i lfc < Ch\\v\\ 2 . 
Proof. Let if> h e P . We have 



(i7 Po (Ag)-A /t (iTp v),V/ l )= E 1*1 (^(A?) - A h (17p v)) 

K£T h 

For all K E Th, using (|2.3[) and the divergence formula, we get 
\K\ (7T Po (Av)-A h (i7p v))| K 

= E ( / Vv ' n ^ CT rfcr _ r,T ( v ( Xi -) - v ( x ^)) ) • 

Thus, by writing the sum over the triangles as a sum over the edges, we get 
(7T Po (Av) - A h (n Po v),^ h ) = fe„-^J^ 



with R ff = J a (Vv ■ n Ki7 . IJ - j- (v(x L J - v(x^J)) da. We denote by D a 

the quadrilatere defined by ~x.k„ , ~*-L a and the endpoints of a. Using a Taylor 
expansion and a density argument, we get as in [7] that 

/ 2 \ V2 

\R°\<Ch[j D ^|H(« 4 )(y)l 2 dy • 
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Thus, using the Cauchy-Schwarz inequality, we get 
|(77p Q (Ag)-A ft (^ Po v),Va)l 



1/2 



< Chi \^L a -^K„ 



EE / |HK)(y)l 2 rfy] 



1/2 



According to ([3]) 



Thus 



the result. 



^7 Po (Av)- A ft (7T Po v),V ft ) < C/i||^|U||v|| 2 . Using ® we get 



5 Convergence of the scheme 

We first recall the stability result that has been proven in [15]. We deduce 
from it an estimate on the Fourier transform of the computed velocity (lemma 
[1]). Using a result on space Po, we infer from it the convergence of the 
scheme (theorem [5]). One shows that if the data uo et f fulfill a compati- 
bility condition [11], there exists a solution (u,p) to equations (P)-© such 
that u € C(0, T; H 2 ) , Vp e C(0, T; L 2 ). We assume from now on that there 
exists C > such that 

(HI) |uO-uoH-i||u^-u(*i)||oo+bi-p(*i)| <Ch, \ul-u° h \ <Ck. 
Let us recall the following result [15] . 

Theorem 1 M^e assume that the initial values of the scheme fulfill (HI). 
There exists a constant C > such that for all m £E {2, . . . , N} 

ki + fcE 11*211* + fc i< - u r _1 i +iEk- ar'iifc < 

n=2 n=2 

and 

m 
n=2 

From now on we set = for the sake of conveniance. One deduces from 
hypothesis (HI) [14] that \p\\ < C and \\u\Wh < C. Now, let e = max(/i, k). 
We study the behaviour of the scheme as e — > 0. We define the applications 
u £ : M -> P , u e : R -> P , u= : R -> P , p e : R -> Pf° and f £ : R -> P as 
follows. For all n G {0, . . . , N — 1} and all £ € [£„, £ n+i] we set 

u £ (£) = u^, u e (£) = < +1 , u°(i) 
p s (t)=p%, f e (£) = 77 Po f (£ n+1 ) , 



= u£ + - (£ - £„) K +1 - fig) , 
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and for all t [0, T] we set u E (t) = u E (t) = u c e {t) = f e (i) = 0, p e (i) = 0. We 
recall that the Fourier transform v of a function v e L X (R) is defined by 

VreR, v(r) = [ e~ 2inTt v(t) dt. (19) 

We have the following result. 

Lemma 1 Let < 7 < j. There exists C > swc/i i/iai /or all e > 

/ |T|^|u £ (r)| 2 dr<C. 

Proof. Let x/ be the characteristic function of an interval I C M. We define 
the application g £ : R — > Po as follows. For all t g" \t\,T\ we set g £ (t) = 0. 
For all t S g £ (i) g P is the solution of 

A ft g £ = A h u e +f e - b h (2u £ -u e (t- k),u s ) + P £ 

with P £ = -V /lPe - 2 fi «(*-fc)-u e + Sifctefcfe) X[t2 T] . W e have omitted 
most of the time dependancies for the sake of concision. Let us estimate g £ . 
We have 

-(A hSe , g £ ) = -(A h u £ , &) - (f £) g £ ) (20) 

+ b h (2 U £ - U £ (t - k) , U £ , g £ ) + (P £ , gfe). 

According to proposition [4] we have 

-(A h g £ ,g £ ) = ||g £ ||£, -(A, 1 u £ ,g £ ) < ||u £ ||h ||g £ ||h. 

Using the Cauchy-Schwarz inequality and |6]) we have 

-%,^)<%\\^\<C\t e \\\ Es \W 
According to proposition [2] and theorem [1] 
|b ft (2u £ - u £ (t - fc),u £ ,g £ )| 

< C|2u £ - u E (t - k)\ \\u E \\ h ||g £ |U < C||u £ |U ||ge||h. 
Using (|18p we have 

P e= (^- + ^X[t 3 ,T^j VhPe + ^X[t 3 ,T]'VhPe(t-k)-^X[t 3 ,T] V h Pe(t-2k). 

Using proposition [3] and the Cauchy-Schwarz inequality we get 

|(P £ ,g £ )| < C (\p £ \ + X [t 3 ,T] \Pe(t - fc)| + X[t 3 ,T] \Pe{t - 2*)|) |div h g £ |. 

Using proposition [5] we have 

|(P £ ,ge)| < C (\ Pe \ + Xlt 3 ,T] \Pe(t - k)\ + X [t 3 ,T] M " 2fc)|) \\z E \\ h . 
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Let us plug these estimates into ([20]) . By simplifying by ||g e ||^ and integrating 
from t = t\ to T we get 



[ T \\Se\\hdt<C + C[ [ T \p s \dt+ C 

Jti \Jt! Jti 



f I 'II -h I \\u s \\ h dt 



< C. 



According to the Cauchy-Schwarz inequality and theorem [T] 

£ \ P e(t)\ dt<VT ^jf \ Pe {t)\ 2 dt\ <Vr(kJ2 |KI 2 ) 

Thanks to the stability of 77p for the L 2 norm we have 

r T N N N 

/ %(t)\dt=kJ2\n Po f(tn)\ <fc^|f(t„)| <^Ell f ll C (0,T;L 2 )^ C - 
J° n=l n=l n=l 

And thanks to the Cauchy-Schwarz inequality and theorem [1] 



N 



[ \\nS)\\hdt<VT f \\u £ (t)\\ldt<CkJ2\\ffX<C. 

Jti Jti n=2 

Thus, since ge (t) -0 forte [0,h}, we get £ \\g £ (t)\\ h = £ \\g £ (t)\\ h dt < C . 
Using definition (JTHJ) we obtain finally 

Vrel, \\Se(r)\\ h <C. (21) 

With this estimate we can now prove the result. Since the function u £ is 
piecewise C 1 on R, and discontinous for t = and t = T, equation (TTT)) reads 

l^- 1 2^ it - k) = AhS " + \ ^ S ° - ^ St) - \ - ** ST+k) 



2 k 



' X[0,ti] ~ o 1 X[T,T+k] 



where So, St t , 5t and 5r+k are Dirac distributions located respectively in 0, 
ti, T and T + k. Let r £ R. Applying the Fourier transform we get 

-2iitt (I -I e- 2i * Tk ) 5 £ (r) = Afcg e (r) + a 



with 



2 2 



'3 u* - fig _ 1 - < - 1 2i7rr \ e- 2 - fc - 1 
2 fc 2 fc }k 
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2.T7 \--le- 2 ™ Tk] - 

2 2 



Taking the scalar product with i u £ (r) we get 

|u £ (r)| 2 =i(A h g e (T),S E (r)) + i(a,fi £ (r) 

Let us bound the right-hand side. According to proposition [4] and ((21 
ifA^T),^^ < ||g £ (r)|U||a £ (r)|U<C||fi E (r)|U. 



On the other hand, using theorem [TJ one checks that a is bounded. Thus, 
according to the Cauchy-Schwarz inequality and J6| 



i[a,u s (T)j < \a\\ue(T)\ <C|a|||u e (r)|U<C||u s (r)|U. 
Hence we have 

VreK, |r||a £ (r)| 2 < C||S £ (r)|U.. 

If r =/= 0, multiplying this estimate by |t| 2 t- 1 , we get |r| 27 |u £ (r)| 2 < 
(7| T |27-i ||u e ( T )||fe. Using the Young inequality and integrating over {r e 
K; |t| > 1} we obtain 

/ |rr|fi £ (r)| 2 dr< / |r| 4 ^ 2 rfr + C [ ||fi £ (r)|| 2 dr. 

J|t|>1 "'|t|>1 >/|t|>1 

For |r| < 1 we have |t| 2 ^ |u £ (t)| 2 < |ti £ (r)| 2 < C||u £ (r)||^ thanks to J6]). 
Thus 

/ \r\^\%(r)\ 2 dr< f |r| 4 ^ 2 dr + C f ||S £ (r)|| 2 dr. 
Jr J\t\>i Jm 

Since 47 — 2 < —1 we have J| T | >:L | -7- 1 4 ^ 2 dr < C. On the other hand, thanks 
to the Parseval theorem and thoremQ] 

/ ||fi £ (r)|| 2 dr< / \\Ur)\\ldt<kf2\\<\\l<C 

Hence the result. ■ 
We introduce the following spaces 

H = {v 6 L 2 ; divv G L 2 et v n\ dn = 0} , V = {v G Hj ; divv = 0}. 
We also set 

2 2 
((u, v)) = ^(Vw„ Vuj) , b(u, v, w) = - ^ (vi, u • VtVi) 



i=l 



i=l 



for all u = (ui,M2) € H 1 , v = («i,i>2) € H 1 , w = {w\,W2) € H 1 . We have 
the following result. 
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Theorem 2 We assume that the initial values of the scheme fulfill hypothesis 
(HI). We also assume that the space step h and the time step k are such 
that h < C k a with a > 1. Then we have u e — ► u in L 2 (0, T; L 2 ) with 

ueC(0,T;H)nL 2 (0,T;V), ^ G £ 2 (0, T; L 2 ). (22) 

at 

We also have u(0) = u and for all ip £ Cq°([0, T]) 

VveV, p^(u,v) + ((u,v))+b(u,u,v)-(f,v)j(if = 0. (23) 

Proof. In what follows, sub-sequences of a sequence (v e ) e> o will still be 
noted (v £ ) e>0 for the sake of convenience. All the limits are for e — * 0. 
According to theorem [1] and hypothesis (HI) we have 

N 

ll u *lk(o 1 r i L") = k W + + *E Kl 2 ^ a 

n=2 

We also deduce from (f6j) , hypothesis (HI) and theorem [1] 

N N 

H^I| 2 L 2 (0,T;L 2 ) = k l^l 2 + k E ^ C + C k E W<Wl ^ C - 

n=2 n=2 

A simple computation shows that there exists C > such that 

II^IIl2(o,t ; L 2 ) - C ^h*(o,T& 2 ) - C - 

Thus the sequences (u e ) e>0 , (u e ) e >o and (Ug) e >o are bounded in L 2 (0,T; L 2 ). 
Therefore there exists u £ L 2 (0,T;L 2 ), u £ L 2 (0, T; L 2 ) and u c £ L 2 (0,T;I/) 
such that, up to a sub-sequence, we have 

u e u , u £ u , — s- u c weakly in L 2 (0, T; L 2 ). 

We claim that the limits u, u, u c are the same. Indeed, let us consider u e — u £ . 
Since uj* - u n h +1 = K ~ < +1 ) + « +1 - < +l ) we have 

N-l N-l 

K - ^II 2 l2(0 , T; l 2 ) < 2fc e k - < +1 i 2 + 2fc E k +1 - < l+1 l 2 - 

n=0 n=0 

According to theorem[T]we have k J2n=o K - K +1 1 2 < C Y,n=o k 3 <Ck 2 . 
Thanks to (|18p we also have 

N-l . N-l N-l 

k £ K+ 1 - u^ 1 ! 2 = ^ 3 £ \^M +1 -pD? <cJ2k 3 <ck 2 . 

n—Q n—1 n—1 

Thus || u E — u e 1 1 i2 (-q t-Ii 2 ) ~ * ^ an< ^ u = ^ ne cnec k s m a simililar way that 
u = u c . Now, using the Fourier transform, we prove the strong convergence 
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of the sequence (u £ ) e>0 in L 2 (0, T; L 2 ). We set v £ = u £ — u. Let M > 0. We 
use the splitting 

/ \v E (r)\ 2 dr= [ |v e (r)| 2 dr+ / |v e (r)| 2 dr = if + jf . 

JR J\t\<M J\t\>M 

Let us estimate jf. Since |v e (Y)| 2 < 2 |u £ (t)| 2 + 2 |u(t)| 2 we have 

Jf <2 / |u e (r)| 2 dr + 2 / |S(r)| 2 dr. 
According to lemma Q] we have 



Thus 




Therefore, for all e > 0, we have — * when M — > oo. We now consider 
4 M . Let r e M. Since u e u in L 2 (0, T; L 2 ), we deduce from definition |[T9]) 

u e (r) — »■ u(r) weakly in L 2 . 

For all i € R we have u e (t) S Po- From definition (fT9|l we infer that u E (r) G 
Po- Now, prolonging u e (r) by outside Q, we deduce from lemma 4 in [7] 
that there exists a constant C > such that 

VryeM 2 , |S £ (t)(-+t,)-5 £ (t)| 2 < ^(^^^(H + C/i). 

Using definition f|19|) . the Cauchy-Schwarz inequality and theorem[IJ we have 

\\%(r)\\l<C [ \\u E (t)\\ldt<CkJ2\\u n h \\l<C. 

J ° n=l 

Thus, using the compactness criterium given by theorem 1 in [7], we get 
u £ (t) —¥ u(t) in L 2 . Thus v e (r) = u e (r) - u(r) — > in L 2 . Therefore for all 
A/ > we have if — > 0. Using the Parseval inequality, and gathering the 
limits for and jf, we get 

/ \v e (r)\ 2 dT= f \v e \ 2 dt = f \u e ~u\ 2 dr^0. 

JR JR JR 

We have proven that u e — > u in L 2 (0, T; L 2 ). 

We now check the properties of u. First, proceeding as in [7], one checks 
easily that u e L 2 (0,T;Hj). Now let q G L 2 (0,T;C£°). According to (O we 
have \7 h (n Pr q) -> Vq in L 2 (0,T;L 2 ). Since u e -> u in L 2 (0,T;L 2 ) we get 

(V h (i7p r g),u e ) -> (Vg,u) = -(g.divu). 
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On the other hand, according to propositions Q] and [3l we have for all e > 

(V h (n Pr q),U e ) = -(ITp«o ? ,diVfcUe) = 0. 

Thus we have J Q T qdivudt = for all g G L 2 (0, T; C§°). Since the space 
is dense in L 2 , we get divu = 0. Hence u 6 L 2 (Q,T; V). Let us now check 

the regularity of — . Using hypothesis (HI), (R)]) and theorem [TJ we have 
at 

du' 



dt 



= l\ul-uW + kf:^<C + Ck±^<C. 



' <iu c 

Thus the sequence ( ) is bounded in L 2 (0,T;h z ). Since u e — > u in 



£ 2 (0,T;L 2 ) with u e L 2 (0,T;H 1 ), proceeding as in, we get 

du': du , , m 
~df Y G (°' T;L ) 

and u G C(0,T;H). 

Let us now prove that u satisfies (f2"3")l . For the sake of simplicity, we omit to 
note some time dependencies. According to (fTT|) we have for all t G 

3 du c 1 du c 1 ~ 

r 1 ^(t-k) A h u £ + b h (2u £ -u e (t- k),u £ ) - f £ 

2 dt 2 dt y ' Re y y h 1 

= ~ + g X[t 3 ,T]J Vfep e + jj X[t 3 ,T] V ft p e (t - fc) - ~ X[t 3 ,T] ^hPe(t - 2 fc). 

Let v G V n (Qf) 2 and V € C°°([0,T]) with i/)(T) = 0. We set y h = 7T Po v. 
Multiplying the former equation by %j) and integrating over \t\ , T] we get 

+ / ijjb h (2u £ - u £ (t - k),u £ ,v h ) dt - ip({ £ ,v h )dt 

Jti Jt! 



T 



X ^(V h p £ ,v h )dt (24) 

With X = -X[t 3 ,T-2k] +|X[t 2 ,t 3 ] -|X[ti,t 2 ] - I X[T-fc,T] - |X[T-2fc,T-M- We 

now check the limits of the terms in this equation. First, according to ([5]), we 
have V/j — > v in L 2 . We will use this limit in the computations below without 
mentioning it. Since if)(T) =0we obtain by integrating by parts 



and 

T ( du c \ r T - k 

^(-^(<-*),v h )dt = -V(0)(ug,v h )- / + A) (u£,v fc )<ft. 
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According to hypothesis (HI) we have u° = u° — * Uo in L 2 and 
u h -> u in L 2 . It implies that (u£,v h ) -> (u ,v) and ^(tj.) (u£, v^) 
^(fc) (u^,v/j) — > ^(0) (uo,v). On the other hand 

t/j'(u c e ,v h )dt = [ X [t u T]^'{^h)dt^ [ i/j'(vL,v)dt 
ti Jo Jo 

and since X[o,T-fc] ^'0 + k) —> ip' in L°°(0, T) 

T-fc 



^'(t+A;)(u|,Vfc)dt= / X[o,r-*]^ / (t+*)(u|,v h )dt-» / ^'(u,v)dt. 
Jo Jo 

Thus we have 

p Gf - 5 f ( *- ^) * - - m (u ' v) - jy Mdt - 

(25) 

Let us now consider the discrete laplacian. Using proposition [6] and the split- 
ting A h v h = (A^Vft - 7T Po (Av)) + 7T Po (Av) we have 

/ i/}(A h u e ,v h )dt = A e + B e 
Jo 

with A e = ffi/, (u e , A, l (7Tp v)-i7 Po (Av)) dt, B e = £ V (u E) n Po (Av)) dt. 
Since 

\A e \ < || A ft (7Tp v)-77 Po (Av) / ^KIM*, 

Jo 

using proposition [7] and the Cauchy-Schwarz inequality, we get 



\A e \<Ch\\v\\ 2 I iP\\u e \\ h dt<Ch I / V 2 ^| 1/ Ml dt 



1/2 / nT , \ 1/2 

e II ft. ' 



1 /2 

Therefore, using theorem [TJ |A e | < Ch (fc£^ =1 WKWh) ^ Ch - Hence 
A £ — > 0. On the other hand, using an integration by parts, we have 



B e = [ ip(u E ,Av)dt-> [ V (u, Av) dt = - [ if) ((u, 
Jo Jo Jo 



v))dt. 



By gathering the limits for A e and B e we get 

I {A h u e , v h )dt -> - / V((u,v))df. 
Jo Jo 

Let us now consider the pressure. We use the splitting 
(VftPe, Vft) = (VftPe, Vft - v) + (VftPe, v - 77 RTo v) + (V 1 h p e , 27rt v). (26) 



16 



Sebastien Zimmermann 



First, integrating by parts, we have 

(VhPe, n RTo v) =-(p E , div (7T RTo v)) + Y] Pe ( i7RT o V ' n K.a)- 

JdK 



KeT h 



Since divv = 0, using the divergence formula and definition (|13jl . one checks 
that div(7TRT v ) = 0. Thus — (p e , div (77rt v )) = 0. On the other hand 

Y / (#aT v-nK,«r) = (( n R'ToV)<T-n Kata ) (p £ \ La -p s \ Krr )da 

KeTh JdK ae£ t nt J a 

and since p e £ Pf c we get Y,KeT h JokP? ( ff RT v • n Ki(r ) = 0. Thus the last 
term in ([26]) vanishes. To bound the other terms, we use the Cauchy-Schwartz 
inequality together with estimates ©, (fT?| and theorem [TJ We get 

|(Vfcp e ,v- v„)| + |(Vhp e ,v- 77rt v)| < Ch\V hPe \ ||v|| 2 < c| ||v|| 2 . 

Plugging these estimates into (PS)) we get J" t T (y h p e , v/j)| di < C ^. By 
hypothesis we have £ < with a — 1 > 0. Thus for e = max(/i, fc) — > 

T 

Xip(VhPe,v h ) dt -> 0. 

1 

Let us now consider the convection term. We set u 6 = 2u £ - u f (t- k) 
and want to find the limit of ip b^(u e , u £ , v/j) dt. We use the splitting 
-b h (u e ,u e , v h ) + b(u, u, v) = A\ + A\ + A| with 

2 

Af = b(u-u e ,u,v) , A% = b(u £ ,u,v) - ^ ( di vK"Q> u l) , 



A% =^(div(w 4 u e ),vl) - b h {u £ ,u e ,v h ). 

i=l 

By definition A\ = — 52i=i ( u i> ( u — u E ) ' Vu$). Using the Cauchy-Schwarz 
inequality we get 

i>\A\\dt< ||V||oo||v|| w i- ||u|| i2(0T;L2) ||u-u £ || i2(0iT . L2) . 



Since u £ — > u in L (0,T;L ) we also have ||u — u e || L2 ^ Q T jjt. — » 0. Thus 
J^tp A\ dt — ► 0. Let us now consider Af. Since u £ • n|ai7 = we obtain by 
integrating by parts b(u e ,u,v) = ( v h div(ui u E )) . Thus 



A 2 = Y ~ v h> div( Ui u £ )) = ^ (vi - v\, u e • Vi 
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Using the Cauchy-Schwarz inequality we get 

f \p\Al\dt < C ||v - v h || L - K|| i2(aT;L2) ||u|| L2(0iT;H i). 

Using a Taylor expansion one checks that ||v — VhJ|i °° < ||v||-w ,1 .° Q h. We 
recall also that ||u £ || L2 , t .j_2n < C. Therefore J^ipA^dt — > 0. Let us now 
bound A%. For all triangle K £ Th and all edge a £ £k H we set 

~ £ ju £ \ K if (u E • TlK,o) > 

K > L " \u E \ L „ if (u £ -n^) <0 ■ 

Using the divergence formula one checks that 

^3 = E E WK ' / ( u ~ J ( n e ' n K,a) da. 
By writing this sum as a sum on the edges we get 

A l = E " v O ' / ( u - *k„,l„) (u £ • n Kattr ) da. 

Thus, using definition (jlip and a quadrature formula 

= J! (v K „ -v L J\a\((n F ncu)(x a ) -u%; aiL J (u £ -n Kat<T ) a . 



|2 



We have \a\ < h and, using a Taylor expansion, one checks that -v^ | < 
h Hvllw 1 •<*>■ Thus, thanks to the Cauchy-Schwarz inequality, we get 

/ \ 1/2 / x 1/2 

\A\\<ChA E |u e (x CT )| 2 E l(^Pru)(x CT )-5' 

Using (HJ) we get 

i^i<c|E^ E m^) 12 

E^ E i(/r Pf .u)(x <r )-3.w a 

v KeT h a££ K n£ z h nt 

Therefore, using a quadrature formula 

|A|| < C\Se\ \n Pr u-u e \ <C\u e \ |iT Pr u-u e | 



1/2 
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By writing TTpn^u — u £ = (ilp^u — u) + (u — u £ ) and using (fT2|) . we get 
I -^-3 1 < C | u E | (h ||u||i + |u — u E |). Thus, using the Cauchy-Schwarz inequality, 
we have 

f M \A%\ dt<C ||u e || L2(0 T . L 2 ) (h ||u|| i2(0iT;H i) + ||u- u £ || L2(0 T . L2) ). 

Thus ipA^dt — > 0. By gathering the limits for Af, A|, A| we obtain 
V' (bh(u E , u e , v ft ) - b(u, u, v)) dt —> 0. Since J* 1 b(u, u, v) dt —> 0, we get 

/ tp (b h (u e ,u £ ,v h ),dt -> / ^>b(u,u, v)di. 
Jti Jo 

Finally, since G Pq, we have: (f e ,Vfc) = (77p f, v^) = (f, v^). Therefore 
i>{f e ,v h )dt= / X[ti,T]V'( f .v/,)d*-» / i/>(f,v)dt. 



We now gather the limits we have obtained into (|24p. The space Vn(C °°) 2 is 
dense in V. Hence we obtain for all v G V and -0 G C°°([0, T]) with ip(T) = 

-V(0) (u , v) - / ^ (u, v) dt + [ V (((u, v)) + b(u, u, v) - (f , v)) dt = 0. 
Jo Jo 

Taking if) = (j> G Cg°([0, T]), we have (f>(0) = and from the definition of the 
derivative in the distributional sense 



pT f^d 

J ^(u J v)* = -y ^-(u,v) 



f/f. 



Thus we have proven (|23|) . At last, let us show that the initial condition 
holds. We have proven before that 

^1^^ weakly in i 2 (0,T;L 2 ). 
dt dt 

Let v G V n (C£°) 2 and i/> G C°°([0, T]) such that ^(T) = 0. We have 
T ( f3 du c F 1 dig. \ , /" T , /du \ , 

Integrating by parts the limit we get 

? T V> (| ^ " ^ (* - fc) , v h ) dt - - V(0) (u(0) , v) - J* i>> (u, v) dt . 

By comparing this limit with (f2"5")) . we get ^(0) (u(0) — Uo,v) = for all 
tp G C°°([0,T]) with ip(T) = 0. Therefore u(0) = u . At last, note that we 
have proven so far the convergence of a sub-sequence of (u e ) e >o towards u. 
But the application u such that (f22|) . (123|) and u(0) = uo hold is unique f[13j. 
p. 254). Thus the whole sequence (u e ) e >o converges towards u. ■ 
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